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Abstract 

We analyze the spectrum of a self-adjoint operator on a Laakso space 
using the projective limit construction originally given by Barlow and 
Evans. We will use the hierarchical cell structure induced by the choice of 
approximating quantum graphs to calculate the spectrum with multiplic- 
ities. We also extend the method for using the hierarchical cell structure 
to more general projective limits beyond Laakso spaces. 

MCS: 34L40 (primary); 34L16; 54B30 

1 Introduction 

Laakso's spaces were first introduced in pTO^ to give examples of spaces that have 
nice analytic properties of any arbitrary dimension greater than one. Among 
these properties is having a sufficient supply of rectifiable curves connecting pairs 
of points to not give positive capacity to single points. In [2], Barlow and Evans 
present a construction of a projective limit space with sufficient conditions that 
Markov processes can be constructed on the limit space and mention that their 
construction should be able to create Laakso's spaces as well. Then in [TS] that 
assertion was proved in detail as well as the construction of a Laplacian operator 
constructed from both the minimal generalized upper gradients provided in [TU] 
and from a Markov process constructed as in [2] . In this paper we work in terms 
of Barlow and Evans construction because their construction as the projective 
limit of quantum graphs (Sections [2] and [3]) makes describing the action of the 
Laplacian simpler (Section [5|). 

Working on the approximating graphs allows us to numerically compute 
eigenvalues and eigenfunctions of the Laplacian on each graph and from the 
projective system of graphs (Section [3]) an eigenfunction on one approximating 
graph can be extended to a corresponding eigenfunction on a more detailed 
approximating graph so that the spectra of the approximating Laplacians are 
building up towards the spectrum of the Laplacian on the limit space. This is 
proved in Section[6l An account of how the numerical calculations are performed 
is given in Sections |4l [5l We begin with a few words on quantum graphs in 
Section [2] and the spectrum of limit 
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2 Quantum Graphs 



A metric graph is a graph which not only has a metric giving the distance 
between vertices but can also give distances between points along the edges. 
Such a graph consist of a set of vertices and a set of edges connecting vertices. 
A pair of vertices can have more than one edge connecting them, and any vertex 
not connected to an edge can be ignored. Each edge is given a length, and the 
distance between any two points is the infimum of the lengths of the paths 
connecting the two points. Such graphs are best visualized as a series of wires 
connected together at the vertices where interesting phenomena can occur not 
just at the vertices but along the "wires" as well. If a Hamiltonian operator 
is defined on the metric graph as well it can then be called a quantum graph, 
the simplest Hamiltonian is the standard Laplacian. A good survey of quantum 
graphs, their properties, and the properties of the Hamiltonians on them is [6] 
followed up in [7]. 

The Hamiltonian need not be an exotic operator, since a quantum graph is 
primarily a collection of intervals operators can be built out of one dimensional 
operators and have their action remain intuitive. A simple operator, and the 
one used in this paper, is the negative second derivative. 

The domain of this operator are those functions on the graph that are twice dif- 
ferentiable on each line segment and such that the directional first derivatives 
at a vertex sum to zero. This vertex condition is known as (Neumann- )Kerchoff 
matching conditions and force the boundary terms to vanish when using inte- 
gration by parts to check the self-adjointness of the operator. If the degree of 
the vertex is one then it is a boundary vertex and the matching condition re- 
duces to a Neumann boundary condition. We will consider this Laplacian on a 
series of quantum graph approximations to Laakso spaces where the Laplacians 
also approximate the Laplacian on the limit space. 

Hamiltonians are self-adjoint operators not merely symmetric ones, see [1] 
for examples in the physics literature. This means that care has to be taken in 
identifying their domains. For instance, second differentiation on the unit inter- 
val is a common Hamiltonian but the domain has to be restricted by boundary 
conditions, classical ones being Dirichlet and Neumann. However on a quan- 
tum graph there are vertices with many line segments meeting. At these points 
there are more choices for boundary conditions, zero (Dirichlet), zero derivative 
(Neumann), one sided derivatives summing to zero (Kirchoff), and so on ^[7]. 
In Section [5] we discuss which vertex conditions are most appropriate to these 
fractals. 



3 Projective Limits 

The Laakso spaces that we will be analyzing Laplacians on are constructed as 
projective limits. The original construction was given in [2], shown to produce 
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Figure 1: Summary of the Projective System, n > m > 



Laakso spaces in [TS], and background for projective limits in the category of 
topological spaces is given in [TT] and [5]. We start with the definitions of a 
projective system and the projective limit space and then move on to giving the 
spectrum of an operator on the projective limit space. 

Definition 3.1. A projective system of topological spaces, (Fi, 0i+i,i)^g, is 
a family of spaces and surjective maps where the Fi are topological spaces and 
4>i+i,i '■ Fi+i — > Fi are continuous. 

Definition 3.2. The projective limit of a projective system, lim^ Fi C Ili^o -^i 

is a topological space. An element {xi} € lim^ Fi if Xi G Fi and (pi^i-iXi = Xi-i 
for all i > and the maps $j : lim^ Fi — > Fj are all continuous surjections. 

Projective limits possess a version of the "universal property." This means 
that if another space L was another topological space satisfying the diagram in 
Figure [T] then all of the maps "factor" through lim^ Fi inducing a continuous 
surjection from L to lim^ Fi. This gives a sense in which the projective limit 
space is minimal because it is the smallest space which has all of the <i>j sur- 
jective. The diagram showing this is Figure [2l where L is the other candidate 
for a limit space and $j the set of surjections belonging to L. By the universal 
property of the projective limit space we have that = o rj for all J > 
and some rj : L ^ lim,_ Fi a continuous surjection. 

We need not only a topology on the limit space, lim^ Fi, but also a measure. 
In [3] sufficient conditions for the existence of a measure on the limit space are 
given. They are: 

1. For A <Z Fi -measurable that 1a dfJ-i = Jf _^_-^ ° 4'i+i,i dfii+i. 

2. The total masses of the measures fii are bounded. 

Given these two conditions there exists a limit measure fi^o with the property 
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This condition forces the sequence of measures fii to be mutually compatible 
enough that the sequence of measures ^*iJ.i{A) := /I^ o $i dfioo, with the 
integral taken over lim^ Fi , converge to the measure ^oo weakly. 

Fix a projective system of quantum graphs, {fnlj^o with measures, Let 
4>„ : lim^ Fi Fn where lim^ Fi is projective limit of F„. It is possible to pull 
back functions on Fn to functions on lim^ i^^ by pre-composing with $„. I.e. 
if / : F„ ^ M then <^*J := / o $„ : hm^ F, ^ M. 

Definition 3.3. Let An be a self-adjoint operator on Fn with domain Dom(An) 
consisting of functions on Fn that are twice differ entiahle on each line segment 
and and also have the property that at each vertex the one-sided first derivatives 
along each edge incident to the vertex sum to zero. These are the Kirchoff 
matching conditions. 

Definition 3.4. A family of self- adjoint operators, (A„, Z3oTO(yl„)) is increas- 
ing if Dom{An) C Dom{An+i) and An+if = Anf for all f G Dom{An) and 
all n > 0. 

Proposition 3.1. The operators An with domains 

Vn = {/o$„|/ e Dom{An)} 

Acting by Anf = {Anf) o for f G I?„ are an increasing family of self-adjoint 
operators and the operator A is the minimal self-adjoint extension of all of the 
An and has domain 

oo 

Dom{A) =\JVn. 

n=0 

Proof. The proof of this statement is rather lengthy and uses Dirichlet forms so 
we refer to Section 7 of [TS] for the proof. □ 
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The self-adjoint operator A is referred to as the projective limit of An- We 
have 

oo oo 

y P„ = y <S>lDom{An) C Dom{A) 

n=0 n=0 

Densely in the graph norm, — ||m||l2 + ||Am||/,2, because of the projective 
limit construction which defines Dom{A) as the minimal domain containing I?„ 
for all n > such that {A, Doni{A)) is self-adjoint. It can be checked that the 
graph norm closure of IJJ^o that domain. To ease the calculations that 

follow we orthogonalize the spaces T>n as follows: 

V'„ = Dom{An) 
V[ = V'q n Pi 

Where is the orthogonal complement of 2?„_i in Ti = L^{L,ii) and 

ffiJ^Lo^'n — U^o -^n- If '^^'^ '^^'^ ^^^'^ ^^'^ spectrum of the Laplacian, —A, when 
operating on each then the spectrum of —A when acting on its entire do- 
main is the union of its spectrum on the P'„ with multiplicities adding. We 
have proved the following theorem. 

Theorem 3.1. Given a projective system of quantum graphs, if (A, Dom(A)) 
is the projective limit of the operators {An, Dom{An)) then 

oo 

a{A) = y a{An\^,). 

n=0 

Where multiplicities of the eigenvalues add. 

Proof. Above we saw that ©J^g^'n = U^o If the union, UJ^o '^{■^n\x>' )j 
did not account for all of the eigenvalues of A then this would imply that that 
there was at least a one dimensional subspace of Dom{A) which is orthogonal 
to every 2?„ which is not true because A is defined as the projective limit of the 
operators, An. Thus a complete set of eigenf unctions for Dom{A) can be found 
in ®^o^'n- A niore detailed version of this proof can be found in □ 

If the geometry of a projective limit space is well understood then it is 
possible to explicitly determine the orthogonalization, !?'„. This is done in 
Section [HI in the case of Laakso spaces. 

4 Laakso Spaces 

Laakso originally described these spaces in pi5^ with emphasis on showing there 
was a large supply of rectifiable curves and a tunable dimension. Later Barlow 
and Evans created a projective limit construction that could be specialized to 
use quantum graphs that would also produce Laakso spaces. This was verified 
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n = l, j = 2 




n=2, j = 2 

Figure 3: Constructions of the simplest case, j — 2. 
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n=2,j=3 



Figure 4: Constructions of a more general case, j — 3. 

in [15] . We consider only a few examples of these spaces to calculate spectra of 
Laplacians, however we take advantage of two ways of visualizing the process. 
One of the factors in determining the dimension of Laakso spaces is how many 
subintervals the new identifications made at each level of approximation make, 
we consider only the cases when this number, j, is fixed from one level to the 
next and work with j = 2, 3, 4, 5, 6, 7. 

The first visualization of the projective limit construction of approximating 
quantum graphs is the account that Barlow and Evans give in . In this account 
the first graph is a unit interval, call it Fq. To construct Fn from Fn-i take two 
copies of Fn-i (if the dimension is larger than two more copies are needed, see 
[15] for details of how many). With those two copies, identify pairs of points, 
i.e. wormholes, from both spaces to create a single space this will be Fn- It 
is important that the locations of the wormholes are chosen so that wormholes 
used for different n do not occur at the same horizontal coordinate. In Figures[3] 
and[4]onc can see the first few steps for when j — 2 and j = 3. This visualization 
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Figure 5: Method of constructing the space when j — 2 



is very convenient for drawing pictures but not for the numerical approximations 
needed to describe a Laplacian and then to calculate its spectrum. 

Another way to visualize the construction of the approximating graphs to 
the Laakso spaces is to think in terms of their cell structure. In general a Laakso 
space's cell structure is not self-similar but in the cases that we consider, where 
at each step we split every horizontal interval into the same number, j, of 
subintervals, the Laakso spaces do have a self-similar cell structure. Again we 
begin with Fq being the unit interval. Now take 2j copies of Fn-i arrange them 
in two horizontal lines of j copies each. Figure [5] shows an example of how this 
is done. To start with we have two copies of _Fo identified together at their 
mid-points, this is shown by the dots, to form Fi. To construct F2 we take four 
copies of Fi with the endpoints label so as to glue the graphs together as one 
would expect from the picture, but the extra step is to say that in the last part 
of the Figure that the points labeled "1" are identified together and the points 
labeled "2" are identified together. The process would them be continued to 
generate F3. This example is when j = 2 but in other cases one just has more 
instances of this process. This may seem like an inelegant mode of thought but 
it leads to simpler labeling schemes for the nodes when building the incidence 
matrices. 

To describe each of the approximating graphs for the numerical computa- 
tions we use incidence matrices. The incidence matrix of a graph is a simple 
representation of how the nodes in space are connected. An entry of indicates 
no connection while a 1 in row 1, column 2, would represent a connection be- 
tween node 1 and node 2. This would also require node 2 to be connected to 
node 1, as the incidence matrix must be symmetric. An entry of 2 means that 
there are two separate edges connecting those two vertices. Many different inci- 
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dence matrices are possible since the labeling of nodes is arbitrary. We assume 
a labeling natural from the second visualization of the construction process in 
the previous paragraph. The incidence matrix according to the vertex labeling 
scheme we use for the n — 1, j — 2 space is 

" 1 " 

10 

110 11 , 

10 
_ 1 _ 

The incidence matrix for the n = I, j = 3 case is 

" 1 " 

1 

1 1 2 

2 1 1 ■ 
1 

1 

The numbering scheme used for these matrices started in the top left node and 
went down columns of nodes, left to right across the space. It is the same scheme 
as shown in Figures [3] and ID These matrices are needed to create the Laplacian 
matrices in the next section. 

From [lOj we know that our space has a Hausdorff dimension Q given by 

1 + ln(2)/ln(j) = Q. 

Since j is an integer larger than one the set of dimensions we deal with are con- 
tained in [1,2]. The more general construction given in [15j allows for different 
values of j at each step. This is necessary to get all possible dimensions possible 
according to [ID]. In these cases we deal with the sequence {ji}^Q listing the 
choice of j at the i'th step in the construction. 



5 Laplacian 

In this section we describe how we numerically approximate the Laplacian in 
Definition 13.31 On each approximating graph an operator acts on a finite di- 
mensional space of functions and so can be given a matrix representation. The 
matrix of the Laplacian is generated by simple manipulations of the incidence 
matrix and is constructed to follow the second difference quotient formula, 

'/ u(x + h) — 2u(x) + u{x — h) 



However, as most nodes are connected to two nodes on each side, we took 
the average of the second derivative along the two paths. To construct the 
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Table 1: Calculated Values of the first 10 Eigenvalues for j = 2, 3, 4, 5, 6, 7 with 
multiplicity, m, and the Expected value, A. 



Laplacian's matrix representation from the incidence matrix, scale each row so 
that the row sum is —2 and then set all diagonal entries to 2. Now the row 
sum is always zero, ff one considers functions depending only on the horizontal 
coordinate, i.e. a function on the unit interval, this Laplacian is coincides with 
the negative second derivative on the unit interval. This representation also 
imposes Neumann boundary conditions at the endpoints by reflecting the value 
at the points just inside the boundary of the space to the points just outside. 
This makes the "first derivative" zero on all boundary points of the space and 
allows the second derivative to be approximated there. The entire matrix is 
then scaled by where h is the distance between nodes or wormholes in the 
space. This distance decreases with increasing n such that the total width of 
the graph is 1 and h is given by h{n) = The algorithmic description is 

desirable here to indicate how the calculations were actually performed. 
The Laplacian matrix for the n = 1, j = 2 space is 
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Figure 6: Multiplicities of eigenvalues when j is 2, 3, 4, and 5 respectively, n = 4 
except for the (a) where n — 6 



The Laplacian matrix for the n = 1, j — 3 space is 



Ml, 3 = 
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The notation is M„j is the Laplacian matrix for the n'th level approximating 
graph to the Laakso space constructing using the value of j. 

This algorithm was implemented in MATLAB. It first produces the incidence 
matrix from recursive patterns emerging from the labeling scheme for the nodes. 
Then the Laplacian matrix is constructed using the process outlined above. 
Since these matrices arc well structured and very sparse, no row has more than 
five entries compared to thousands of zeroes, we were able to find the thousand 
smallest eigenvalues numerically using the MATLAB eigs command which is 
based on ARPACK, a package that implements an Implicitly Restarted Arnoldi 
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(a) A = 9.74 



(b) A = 37.49 



Figure 7: Eigenfunctions on the n = 3, j — 2 space 

Method. More information on ARPACK can be found in its users guide, [T^]. 
Even using sparse storage techniques the computations are impractical on a 
personal computer much past n — 9 when j ~ 2. 

Fix the Laakso space, L, where j is the number of subintervals in the con- 
struction. The action of the operator Mn(l)njy^,n > m, on / e C{Fm) needs 
to be related to the action of Am to justify approximating A by Mn instead 
of An. First (pn^m ■ -> Frn is & surjcction, so 0* ,„ : C(i^„) C{Fn) by 
'Ki mf — f ° 4'n,m- Sincc M„ is a second finite difference operator on the graph 
of Fn it is seen that Mn4>n m for n > m are the standard second finite difference 
operators on each line segment of F,„ which converge to second differentiation 
as n goes to infinity. Then this holds for functions on any Fm so the Mn can 
be taken as suitable approximations to second differentiation. The operator 
An (Definition 13. 3p acts by second differentiation on each line segment of the 
quantum graph Fm- Then for fixed j and m and / G Dom{Am) 

lim Af„0j; „/ = AnJ. 

n — ^CJO 

The proof of this statement is standard and omitted, but is essentially just a 
statement about approximating derivatives finite difference operators. 

Using the matrix representations of Mn given earlier in this section nu- 
merical estimates of the eigenfunctions and eigenvalues were computed using 
MATLAB's eigs function which calls the ARPACK package, [12]. The ten low- 
est eigenvalues for each of the spaces with j = 2, 3, 4, 5, 6, 7 are listed in Table 
[1] along with their observed multiplicities, and histograms given in Figure [6l 
In the next section we give a theoretical account of how these spectrums with 
multiplicities arise and the calculations to show the exact details for j — 2,3 
as the last result of this paper. Figure [7] show two eigenfunctions labeled with 
eigenvalues for the Laakso space with j = 2. These two pictures show how 
high multiplicity can be obtained since the eigenfunctions are only obeying the 
Kirchoff matching condition at the vertices of the graph. These matching con- 
ditions are weak enough so that functions which do not appear smooth at the 
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vertices are still in the domain of the Laplacian. 



6 Spectrum of The Laplacian 

In this section we use the framework from Theorem 13.11 to verify that the cal- 
culations in the previous section are valid. Fix a Laakso space Lj = L, and 
the sequence of quantum graphs of which it is the projective limit, {Fn}^^Q, 
[21 [TJ] . Let <&„ : L — > F„ where L is a Laakso space and F„ one of the ap- 
proximating quantum graphs be a projection map. It is possible to pull back 
function on Fn to functions on L by pre-composing with I.e. if / : — > M 
then / := / o $„ : L ^ M. We can then take A — Aj to be the self-adjoint 
operator with it's domain Dom{A) where this operator is the projective limit of 
An which acts by second differentiation on each line segment of F„ with domain 
Dom{An) as given in Definition 13.31 Recall also the orthogonalization of the 
domains of An as 2?^ from Theorem 13. II 

In the case of the Laakso spaces it will be useful to have a clear picture 
of what the projection maps </>„ do to a function from C(i^„_i) as it is being 
pulled back to an element of C(F„) to better understand what the orthogonal 
space is. Since F„ is constructed from copies of F„_i modulo identifications, 
we can think of Fn = Fn-i x G„ so a function from C(-F„_i) when pulled back 
to C{Fn) will be constant across Gn for a given coordinate in Fn-i- However, 
at wormholes G„ is identified to a single point so the pull back of a C(F„_i) 
function will just have its value there. If we want to describe V'n — T)'^^_^f\'Dn- 
we have to describe the functions on Fn that are orthogonal to those on Fn-i- 
The functions pulled back from C{Fn-i) to C{Fn) are constant across G„ so the 
orthogonal functions would be those that average to zero across G„ and at the 
n — th level wormholes where G„ is identified to a single point the orthogonal 
functions must be equal to zero. 

In [TU] a sequence of integers {ji} is associated to each Laakso space in the 
examples we calculated we have taken j — ji but we use the Laakso's notation 
for the generality since nowhere in the following proof do we use the fact that 
the ji are equal. Let A — Aj be the self-adjoint operator for the Laakso space 
with ji = j, and An = Anj the self-adjoint operators on the approximating 
quantum graphs. 

Theorem 6.1. Let the sequence dn = YVj^iJ^^ n > be associated to a given 
Laakso space, with Hausdorff dimension, Q, between one and two. Set do — I. 
Then the spectrum of A on this Laakso space is 



^ f , 2_2 ^ oo oo . ,22 ^ oo oo 

-(^)=UU{^}uUU{^}uUU 

n=Ofc=l n ■> n=2k=l ^ " ^ n=lk=0 



(2fc-H)27r^ 



Proof. The proof falls naturally into two parts. The first is to show that An 
restricted to V'n has the claimed spectrum. The second is to show that the 
union of cr(A„) contains all of a (A). Calculating exact multiplicities is left for 
the next theorem. 
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Figure 8: The three types of pieces that the orthogonality condition creates in 



The approximating graphs for Laakso spaces are introduced in [15j and have 
been reviewed Section 3) For example if ji — 3 then we have the graph in 
Figure 2] as the second approximating quantum graph. The structure that we 
exploit to calculate the spectrum is that the graph can be broken down into 
simple configurations of line segments on which the behavior of An_j is the 
usual second differentiation. All of the operators, An, have Kirchoff matching 
conditions which imposes Neumann boundary conditions since the boundary 
consists of degree one vertices. 

The functions in 2?'„ are those functions in the domain of the operator A 
that are expressible as $* / where / G Dom{An) where the values of / summed 
over the copies of F„_i that F„ is built from are all zero, especially at the n — th 
level wormholes where all the copies of have been collapsed to single points 
/ must be equal to zero at these wormholes. Referring back to Figure SI if one 
considers the two most upper left line segments forming the sideways "V" . An 
element of 'D'2 would have to equal zero at the second level wormhole, i.e. the 
apex of the V. It must also sum to zero over the copies of Fi from which this 
graph is made meaning that the value of / on the upper branch of the V must 
be opposite to the value of / on the lower branch. On this "V," the part of 
the function that is symmetric across the two branches is in 2?'i, while the 
anti-symmetric part is in 2?' 2 CT>2- 

The zeros at the n'th level wormholes break the graph F„ into three kinds 
of pieces. The first are those segments at the boundary which have length cJ„, 
Neumann boundary conditions on the boundary end and Dirichlet boundary 
conditions on the interior end. If jVt > 2 then we have interior intervals, such as 
those that form the loops in Figure 31 which again have length dn butthese have 
Dirichlet boundary conditions at each end. The third configuration is a cross, 
four segments of length dn joined at a common point and Dirichlet boundary 
conditions at the four outer vertices. These three configurations can be seen 
in Figure [H In the cross we must note that the upper and lower branches are 
from different copies of Fn-i so / is not forced to equal —g for functions in 2?'„ 
but we will consider the symmetric and anti-symmetric cases where f ~ g and 
/ — — g from which any function on the cross is a combination, li f ~ g then 
the cross is essentially just an interval of length 2dn with Dirichlet boundary 
condition and if / = —g then at the intersection / = and we have two intervals 
of length dn with Dirichlet boundary conditions. 



Fn. 
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Note that the crosses do not appear for n = or 71 = 1 for any jo or 
ji. The intervals with Dirichlet conditions do not appear in n — or when 
jn ~ 2. The intervals with mixed boundary conditions do not appear when 
n = either. What does appear when n = is an interval of unit length 
and both ends having Neumann boundary conditions. When all of these are 
put together and their spectra found by the usual methods we get the claimed 
spectra where the contributions towards multiplicity from each of the pieces 
may be zero occasionally, i.e. ji = 2 so there are no intervals of length dn and 
Dirichlet boundaries except those which are part of the crosses. 

The last claim that these are all of the eigenvalues of A on L follows from 
the fact that IJ^^Lo -^''i dense in Dom{Aj) by the projective limit construction 
of the operator A. So if there were another eigenvalue unaccounted for there 
would be at least a one dimensional subspace of Dom{A) orthogonal to a dense 
subset of Dom{A), since this can't happen the entire spectrum is now accounted 
for. □ 

There is a more general statement for higher dimension that replaces the 
crosses with 2*^+^ line segments of length c?„ on each side of a central vertex 
where k is the number of Cantor sets in the construction of the Laakso space 
as in [10] . The spectrum given in these higher dimensions is given by the same 
expression but the multiplicity counting as in the following Proposition has to 
be redone for each k. 

Proposition 6.1. The first ten eigenvalues of Aj and their multiplicity for the 
j = 2 and j = 3 constructions determined from Theorem \6.1\ agree with the 
computed values in Table{^ 

Proof. It is immediate from the spectrum of An.j on a given Laakso space as 
determined in Theorem 16.11 that the lowest eigenvalue is increasing with n, so 
to find the first ten eigenvalues only a small number of approximating graphs 
are necessary. 

Let us establish some notation. Let a'^[0,d]^ be the spectrum of the usual 
second derivative on an interval of length d. If either, or both, of the iV's 
are replaced by then Dirichlet boundary conditions are indicated instead. 
Denote by (To[X, djg the spectrum of second differentiation on a cross where 
the four intervals have length d and Dirichlet boundary conditions at the four 
ends points. If we write 2 x (To[X, (i]Q the two denotes the multiplicity of these 
eigenvalues. 

For j = 2, we get from drawing the first few approximating graphs and 
forming the function spaces V'n and calculating the spectrum of Anj on each 
of these spaces we get the following table. Let cr„ be the spectrum of the An j 
operator on the function space 



ai = 2xa^[0,i] 
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02 = 4xa^[0,i]" + lxaS[X,l]!] 
a, = 8xa^[0,i]° + 6xa!][X,i]!] 
a4 = 16xa^[0,^]" + 28xaO[X,^]0 

It is easy to check that 1 x crg[X, d]^ = 1 x^ [0, 2d]^ + 2 x crO[0, d]". Using the 
known spectra for second differentiation on intervals the first ten eigenvalues are 
readily computed along with their multiplicities which do match those found in 
Table [TJ In the case where j = 3 we need fewer approximating graphs to 
account for the first ten eigenvalues. The spectra of the Laplacian on the first 
few functions spaces are: 



do = lxa^[0,l]^ 

<Ji = 2xa^[0,if + lxaO[0,i]° 

a2 = 4xa7V[0,i]O + 3xaO[0,i]O + 2xaS[X,i]O 

When these spectra are written out with the indicated multiplicities and com- 
piled into a single list we get the same eigenvalues and multiplicities listed in 
Table [H □ 
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